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Quantum depletion from an atomic quasi one dimensional 
Bose-Einstein condensate with a dark soliton is studied in a 
framework of the Bogoliubov theory. Depletion is dominated 
by an anomalous mode localized in a notch of the condensate 
wave function. Depletion in the anomalous mode requires dif- 
ferent treatment than depletion without anomalous modes. 
In particular, quantum depletion in the Bogoliubov vacuum 
of the anomalous mode is experimentally irrelevant. A dark 
soliton is initially prepared in a state with minimal depletion 
which is not a stationary state of the Bogoliubov theory. The 
notch fills up with incoherent atoms depleted from the con- 
densate. For realistic parameters the filling time can be as 
short as 10 ms. 

03.75.Fi,03.75.-b,05.30.Jp 



I. INTRODUCTION 

Since first experimental realizations of atomic Bose- 
Einstein condensates (BEC) [1] the condensates have 
been subject to intensive experimental and theoretical 
research. The condensates are relatively easy to manip- 
ulate and a rich variety of phenomena can be studied in 
a well controlled way. One example is a dark soliton in 
analogy to nonlinear optics [2-4]. The dark soliton is 
a place in a quasi one dimensional condensate where a 
condensate wave function <j>(x) changes its sign. Its wave 
function has a notch where density of the condensate is 
zero. This notch appears as a dark spot on a bright cigar 
shaped image of a condensate. Dark solitons were re- 
alized in the experiment of the Hannover group [5] and 
their 3 dimensional analogs were realized in Ref. [6] . 

A wave function of the dark soliton is proportional to 



tanh[(a; — q)/l] 



(1) 



where I is a local value of the healing length at the posi- 
tion q of the soliton. The density of coherent atoms in the 
notch at x = q is zero. This creates a hole that might be 
filled up with incoherent atoms. Given that interactions 
between atoms are repulsive, depletion of atoms from the 
condensate to the notch is energetically favourable. 

Condensate with a soliton at q = has an antisymmet- 
ric wavefunction of x. This means that all N atoms are 
in the same antisymmetric single particle state. The an- 
tisymmetry of the condensate wavefunction is preserved 
by evolution according to the Gross-Pitaevskii equation. 



However, the full iV-body Hamiltonian depletes pairs of 
atoms from the antisymmetric condensate wavefunction 
to symmetric modes of the trap. This two atom pro- 
cess does not change the symmetry of the iV-body state. 
The symmetric modes have nonzero density in the soliton 
notch at x = 0. Simply because the condensate density 
is zero at x — 0, the x — is where we can clearly 
see the incoherent atoms. In our previous paper [7] we 
studied a toy model of the dark soliton. We truncated 
single particle Hilbert space to 3 lowest modes of a one 
dimensional harmonic trap. This is a self-consistent ap- 
proximation for small N. We have shown that a state 
where all atoms are condensed in the same antisymmet- 
ric wave function is not a stationary state of a quantum 
Hamiltonian. Atoms are gradually depleted from the an- 
tisymmetric condensate wave function and fill its notch. 
In this paper we reach similar conclusions in the opposite 
regime of very large N which is relevant for the Hannover 
experiment [5]. We use as a tool the Bogoliubov theory 
with noninteracting quasiparticles. 

In the Bogoliubov theory one considers small fluctu- 
ations around a stationary condensate background. In 
a uniform condensate the soliton (1) has a translational 
mode 9 g tanh[(x — q)/l] with zero frequency. The same 
soliton in an nonuniform condensate in a harmonic trap 
has an anomalous mode with negative frequency equal 
to —l/y/2 times the trap frequency ui x [4]. The nonuni- 
formity breaks the translational invariance, moving the 
soliton away from the maximal density in the center of 
the trap lowers its energy. We will see that the anoma- 
lous mode is responsible for depletion of atoms from a 
condensate to the soliton notch. Depletion in the Bogoli- 
ubov vacuum has been recently studied in Ref. [8] . How- 
ever, we argue that depletion in the Bogoliubov vacuum 
analyzed in Ref. [8] is not experimentally relevant. We 
construct a quantum state with minimal depletion which 
approximates the state of the system right after creation 
of the dark soliton by phase imprinting [5] . Time evolu- 
tion of such a state is responsible for greying of the dark 
soliton. 

In Ref. [3] it is suggested that deviations from the mean 
field description of the dark soliton observed experimen- 
tally [5] can be explained by dissipation due to collisions 
of the condensate with a thermal cloud of noncondensed 
atoms. There is no doubt that dissipation influences dy- 
namics of the soliton in that experiment. However, we 
show that the dark soliton is greying even in the absence 
of any thermal cloud or even in ID where the cloud de- 
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couples from the soliton [3] . The notch fills up with inco- 
herent atoms quantum depleted from a condensate. Our 
calculations, which are idealized with respect to the ex- 
periment, give a greying time of O(10ms), which is con- 
sistent with the experimental observation. 

This paper is organized as follows. Section II gives 
several definitions and approximate formulas valid in the 
Thomas-Fermi regime of strong interactions. Section III 
gives Thomas- Fermi approximate analytic expressions for 
the frequency and the cigenfunctions of the anomalous 
mode. In Section IV we work out a general expression 
for a time-dependent number of atoms depleted in the 
anomalous mode. In Section V the number of depleted 
atoms is calculated in the Bogoliubov vacuum of the 
anomalous mode. In Section VI we construct a quan- 
tum state with minimal depiction which approximates 
the state of the system right after creation of the dark 
soliton. In Section VII we follow time evolution of the 
number and density of depleted atoms, and estimate the 
greying time after which the notch fills up with incoherent 
atoms. Section VIII contains discussion and conclusions. 



II. THE DARK SOLITON 

In dimcnsionless oscillator units a Gross-Pitaevskii 
equation [9] for the condensate wave function <j>(t, x) in a 
harmonic trap is 



id v 



4* 



1 2 



(2) 



<p(t,x) is normalized to 1. The ground state of this 
equation is a stationary state <fio{x) exp(— inot) with a 
symmetric <po, <f>o{%) — (f>o{—x). In the Thomas- Fermi 
regime, g 3> 1 , (f>o can be approximated by 



0o (&) = 



' 2/j, — x 2 



2ff 



(3) 



This wave function is normalized to 1 provided that fi = 
(3g/2) 2 / 3 /2. It is localized and varies on a lengthscale of 



R = 



1/3 



(4) 



the so called Thomas-Fermi radius. 

The lowest antisymmetric stationary state of the 
Gross-Pitaevskii equation (f>i(x) exp(-ifiit) is a dark soli- 
ton localized in the center of the trap. In the Thomas- 
Fermi regime <j>i can be approximated by the ansatz 



Here 



M x ) 



lo = 



S>1 



<Po(x) tanh(x/Zo) 



V2 



(I) 



1/3 



(5) 



(6) 



is a healing length of the condensate at the condensate 
peak density in the center of the trap. The soliton solu- 
tion in Eq.(5) has a notch at x = 0. We note that with 
increasing g the Thomas- Fermi radius R is growing while 
the healing length l tends to zero. For g 3> 1 the width 
of the soliton is very small as compared to the size of the 
condensate, lo -C R. 



III. THE ANOMALOUS MODE 

Bogoliubov equations [9,10] for small perturbations on 
the dark soliton background are 

1 '/ 1 9 2 2 

— -u + —x u + 2g<p 1 u — fiiu + gepiv = +luu , 
~v" + ^x 2 v + 2g(j>\v - mv + g<t>\u = -wv . (7) 

Here = d/dx. We introduce two functions f±(x) = 
u(x) ± v(x) and combine Eqs.(7) into 



(8) 



~\f- + \* 2 f- + 9<frlf- - Mi/- = w/+ • (9) 

Eqs.(8,9) have a solution with negative lo [3,4], the so 
called anomalous mode. In the Thomas-Fermi regime lu 
can be approximated by [3,4] 



LJ ^ 



-1 

71 



(10) 



i.e. —\j\pl times the frequency of the harmonic trap. 
This Thomas-Fermi value of to follows from effective 
equations of motion for a position of a dark soliton in 
Ref. [4], and from a perturbative (g 3> 1) treatment of 
the Bogoliubov equations [3]. 

The anomalous mode is a perturbative approximation 
to the motion of the dark soliton with respect to the back- 
ground condensate. In the limit of g — ► oo the Thomas- 
Fermi radius of the condensate (4) is growing to infinity 
and at the same time the healing length (6) is shrink- 
ing to 0. In this limit the condensate is uniform on the 
lengthscale of the soliton width and the anomalous mode 
of Eqs.(8,9) becomes 



lim f+(x) oc 



1 



cosh (x/lo) 



lim = . 



(11) 



This f+(x) is simply a translational mode of the soliton 
in an almost uniform condensate. This limit gives us the 
leading g 1 term in the expression for / + (x) . 

We use Eq.(9) to find the leading term in f-(x). To 
begin with we note that 4>q{x) satisfies 
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'o + ^% + 50o - Mo^o = 



(12) 



Comparing Eq.(12) with Eq.(9) one may conclude that 
f-(x) = a<j)o(x) with a certain constant a. Substitution 
of this ansatz and of leading f+(x) from Eqs.(ll) into 
Eq.(9) and then using Eq.(12) gives an equation 

-(Mi - Ho)a4>o{x) + 
5 [<f*\{x) - 4>l{x)] a(j) (x) 



cosh (x/l ) 

We expect that with a right choice of a this equation is 
approximately satisfied for g ^> 1. As /ii — /io = O(g ) 
the first term on the left hand side is small as compared to 
the second term and can be neglected. The second term 
is localized on a healing length around x = and the 
slow factor 4>q(x) in this term can be replaced by 0o(O). 
Finally we can use the g » 1 formulas in Eqs.(3,5) to get 



g [tanh 2 (x/l ) - l] a$(0) 
This equation is satisfied by 



cosh (x/Iq) 



5^(0) 3^/5 



(13) 



(14) 



which defines the leading term in /_ (x) . After normal- 
ization, so that / dx (\u\ 2 — \v\ 2 ) = J dx /+/_ = 1, we 
get 



/+(*) 
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2\/2 cosh 2 (x/l ) 



9>1 2 

f-( x ) - \J^M X )- 



(15) 
(16) 



We note that f+(x) is localized within the notch of the 
dark soliton (5). 



IV. NUMBER OF DEPLETED ATOMS 

We will argue that depletion of atoms from the conden- 
sate fills up the notch in the condensate wave function 
(5) with incoherent atoms. This effect will appear exper- 
imentally as greying of the dark soliton. Solution of the 
equations (7) reveals many different modes corresponding 
to the Bogoliubov spectrum. However, according to the 
numerical results in Ref. [8], only the anomalous mode 
has a wave function localized in the soliton notch [8], 
compare Eq.(15). From all the modes the anomalous 
mode will contribute the most to the density of incoher- 
ent atoms filling the notch. This expectation is confirmed 
by numerical calculations of incoherent atoms density 
distribution in the Bogoliubov vacuum state [8], where 
the density, which is strongly peaked in the notch, is 
dominated by a contribution from the anomalous mode. 



In the following we truncate the Bogoliubov spectrum to 
the anomalous mode alone. 

Having constructed f± = u ± v we can calculate an 
operator corresponding to a number of depleted atoms 
in the anomalous mode 

dN(t) = J dx 5ft (x,t) 84>(x,t) , (17) 



where 



6<f>(x,t) = u{x) e~ luJt a + v*(x) e +iwt a f 



(18) 



is an annihilation operator of an incoherent atom in the 
anomalous mode. The operators a, annihilate/create 
a Bogoliubov quasiparticle and they fulfill the usual 
bosonic commutation relation, [a, a'] — 1. We use 
Heisenberg picture defined by the Bogoliubov Hamilto- 
nian of the anomalous mode 

H B = -\u\a)a . (19) 

The operator of the number of depleted atoms is 

dN(t) = J dx [\u\ 2 a)d+ \v\ 2 da) + 

uve- 2lut aa + u*v*e +2wt a)cJ] . (20) 

The density of incoherent atoms is the expectation value 
of the integrand. 



V. DEPLETION IN THE BOGOLIUBOV 
VACUUM STATE 

For a condensate in a ground state like in Eq.(3) the 
Bogoliubov spectrum has only positive frequencies (i.e. 
there are no anomalous modes) . The state with no quasi- 
particles, the Bogoliubov vacuum, is the ground state 
of the Bogoliubov Hamiltonian. In the absence of any 
anomalous modes it is natural to calculate the depletion 
in the Bogoliubov vacuum - the T — state of thermo- 
dynamic equilibrium [10]. 

Following the numerical calculation in Ref. [8] we cal- 
culate analytically for g 3> 1 a number of depleted atoms 
in the vacuum state of the anomalous mode i.e. the state 
|0) a annihilated by a, a\0) a — 0. This vacuum state is 
an eigenstate of the Bogoliubov Hamiltonian (19) so the 
number is stationary 

dN a = a (0\dN(t)\0) a = j dx \v{x)\ 2 , (21) 

and scales with g as dN a oc g 2 ^ 3 . We compare the density 
of depleted atoms in the soliton notch 



dNa 

lo 



oc g oc N , 



(22) 
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(where N is the total number of atoms in the system) to 
the density of atoms in the condensate 



N N ,t2/3 



R 



,1/3 



(23) 



We see that, for a sufficiently large total number of atoms 
N the density of incoherent atoms in the notch will ex- 
ceed the density of the background condensate. The 
notch will not be visible. When the density of depleted 
atoms is comparable to the density of the condensate one 
cannot neglect interactions between quasiparticles. How- 
ever, because of the potential instability which shows up 
in the negative frequency of the anomalous mode, we 
qualitatively expect that, unlike in more common situa- 
tions where all Bogolubov frequencies are positive, colli- 
sions between quasiparticles will not suppress depletion. 

In Fig.l we show condensate and total density plots 
in the Bogolubov vacuum for the parameters of the Han- 
nover experiment [5] (N — 1.5T0 5 87 Rb atoms in the har- 
monic trap with uj x — 2ir ■ 14 s _1 and lo± = 2ir ■ 425 s _1 ; 
g = 7500 [11]). As we can see, even for very reason- 
able trap parameters and moderate N the soliton notch 
is substantially filled with incoherent atoms. For a larger 
N, the soliton would not be observed at all. 
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FIG. 1. Single particle density function corresponding to 
the experimental parameters [5]. The coherent part |0i(a:)| 
[Eq. (5)] is depicted in panel (a), the incoherent part |w(a:)| 2 
corresponding to the vacuum of the anomalous mode |0) a is 
shown in panel (b). Panel (c) shows the total (i.e. the sum 
of the coherent and incoherent parts) density function and 
panel (d) is a zoom on the x ~ part of panel (c) — note that 
the density of incoherent atoms at the soliton notch is about 
14% of the maximal density of the surrounding condensate. 

The depletion in the Bogoliubov vacuum of the anoma- 
lous mode is not relevant for dark solitons in current ex- 
periments 

• The vacuum |0) a is not a ground state of the Bogoli- 
ubov Hamiltonian (19). Creation of a quasiparticle 
in this vacuum with the creation operator actu- 
ally decreases the energy. The Bogoliubov vacuum 



is not the state of thermodynamic equilibrium at 
T = 0. 

For large enough TV the density of incoherent atoms 
in the state |0) a exceeds the density of the con- 
densate. The linearized Bogoliubov theory breaks 
down and the calculation of depletion in the state 
\0) a is not self-consistent. 

The Bogoliubov vacuum |0) a is not the state with 
minimal possible depletion. When v(x) ^ 
quasiparticles cannot be identified with incoher- 
ent atoms, compare Eq.(18), and the state with 
no quasiparticles |0) a is not the state with mini- 
mal number of incoherent atoms. It is possible to 
construct a state with a very small depletion as 
compared to the depletion in the state |0) a , see the 
next Section. 



VI. THE MINIMAL DEPLETION STATE 

Solitons are excited experimentally employing a phase 
imprinting method [5] . We assume that before the phase 
imprinting the system of N atoms is in its ground state. 
The ground state is a condensate (3) with all Bogoliubov 
modes in a vacuum state. There are no anomalous modes 
and the density of incoherent atoms is very small as com- 
pared to the density of the condensate. The ground state 
is very close to a perfect condensate with all atoms in 
the condensate wave function (3). We approximate the 
initial condensate by a perfect condensate without any 
depletion. 

For a perfect initial condensate the phase imprinting 
can be thought of as a nearly instantaneous operation 
on the condensate wave function. Right after this oper- 
ation one gets a condensate with an antisymmetric wave 
function <f>(x) = —(j){—x). To focus on the quantum de- 
pletion and not on classical evolution of (f>(x) we assume 
that 4>{x) = (f>i(x). The initial dark soliton is a (nearly) 
perfect condensate with the wave function (5). In the 
Hannover experiment [5] they create moving grey soli- 
tons while here, for the sake of simplicity, we assume a 
stationary dark soliton. We think that in spite of this, 
qualitative comparison with the experiment is still pos- 
sible. 

The Bogolubov vaccum |0) a is an eigenstate of the N- 
body Hamiltonian. For large N the |0) a has substantial 
depletion. The initial state without depletion is far from 
|0) a or from any other eigenstate. The initial state is not 
stationary, its depletion may grow from zero to values 
which far exceed the depletion in the Bogolubov vacuum. 
In the following we quantify this effect. 

In the framework of the Bogoliubov theory truncated 
to the anomalous mode (19), the initial depletion-free 
soliton can be very well approximated by a variational 
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state IV') that minimizes the initial number of atoms de- 
pleted from the background condensate (j)\{x) 



dN(0) = (ip\dN(0)\ip) . 



(24) 



To find this optimal state it is convenient to introduce 
bosonic operators b and b < , 





where 



(a + a f ) 



n = 2-^i(3 5 ) 1 / 3 



(a -at) , (25) 



(26) 



In terms of these new bosonic operators the initial num- 
ber of depleted atoms (20) becomes 



This operator is a sum of a nonnegative operator ^ wb 
and a constant, wb > is minimized by the unique vac- 
uum |0)b of the bosonic annihilation operator b, b\0) b = 0. 
The minimal depletion state is \ip) = |0){,. This vacuum 
is different from the Bogoliubov vacuum |0) a , 



|0>6 ^ |0)„ . 



(28) 



In this optimal state |0){, the initial number of depleted 
atoms in the soliton notch scales with the number of 
atoms as TV 1 / 3 . The initial density of incoherent atoms 
in the notch 



^(°) „ ^2/3 



oc 



(29) 



where dN b {0) = 6 (0|aW(0)|0) 6 , grows with N in the 
same way as the condensate density in Eq. (23). 

In Fig. 2 (a,b) we show that for a realistic value of the 
scattering length and other parameters the minimal den- 
sity of incoherent atoms in the notch is around 1% of the 
peak density of the surrounding condensate. Our mini- 
mal depletion state is a very accurate approximation to 
the initial perfect condensate with a dark soliton. As 
both coherent and incoherent densities grow with N in 
the same way compare Eqs. (23,29), the quality of this 
approximation does not change with increasing N. The 
negligible depletion in the minimal depletion state jus- 
tifies a posteriori our truncation of the theory to the 
anomalous mode. 

Moreover, close to the initial minimal depletion state 
|0)b we can safely use the noninteracting Bogoliubov the- 
ory. This theory assumes that the density of incoherent 
atoms is small as compared to the density of the con- 
densate. In contrast, for large N the Bogoliubov vacuum 
\0) a does not satisfy this assumption. 
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FIG. 2. Panel (a): the initial (t = 0) single particle den- 
sity function for the system initially in the minimal depiction 
state 1 0)b (the values of the parameters correspond to the 
Hannover experiment [5]). Panel (b) is a zoom on the x ~ 
part of panel (a) — the density of the incoherent atoms at the 
soliton notch is about 1% of the maximal condensate density. 
Panels (c)-(d): the same as in the previous panels but for 
t G = 0.42 (tt/\/2wx) = 10.6 ms. 



VII. TIME EVOLUTION OF THE NUMBER OF 
INCOHERENT ATOMS 

The initial density of incoherent atoms in the notch 
(calculated in the |0)h state) is negligible, see Fig. 2. The 
initial dark soliton is very dark. This small depletion is 
not stationary because the initial minimal depletion state 
1 0) (, is not an eigenstate of the Bogoliubov Hamiltonian 
(19). From Eqs. (20,25) one can easily deduce time evo- 
lution of the number of depleted atoms 



dN b (t) 



6<0|dW(t)|0> 6 

(fi 2 - 



dN b (0) + 



- 2 ) 



2\2 



sin 2 (uit) 



(30) 



The coefficient in front of the sin 2 (cut) function scales 
with the total number of atoms as ./V 1 which implies 
that the maximal density of depleted atoms in the soliton 
notch 



dN, 



(*- zr) 
lo 



AT4/3 



(31) 



Hence, for large enough TV the density of depleted atoms 
in the soliton notch will reach the density of the con- 
densate (23) before it reaches the maximum (31) at 
t = tt/2|w|. 

Eq. (30) might suggest that the system will go through 
a series of periodic revivals of the initial state \0}b with a 
period of 7r/|w|. This is not a correct conclusion for large 
N. After a greying time 



t G < 



2U 



(32) 
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when the density of incoherent atoms in the notch and 
the density of the condensate become comparable for the 
first time, we cannot use the noninteracting Bogoliubov 
theory any more [10]. After to one can ignore neither 
anharmonicity in the Hamiltonian (19) of the anomalous 
mode nor interaction of the anomalous mode with other 
Bogoliubov modes (phonons). We expect that these non- 
linearities lead to fast dephasing and prevent any revivals. 
The number of depleted atoms in the notch initially fol- 
lows Eq.(30) until at around the greying time to their 
density saturates at a density comparable to the density 
of the condensate. This will appear as greying of the 
dark soliton. The nonlinear stage deserves more quanti- 
tative analysis. However, we expect that, because of the 
global instability that shows in the negative frequency of 
the anomalous mode, the nonlinearities will not suppress 
depletion. 

Eq.(32) gives a crude upper estimate of to, to < 
n/2\u)\. Given Eq.(10) we obtain in real time units 



to < 



V^UJx 



(33) 



Here lv x is an axial trap frequency. In the Hannover ex- 
periment of Ref. [5], w/\/2u x — 25 ms while the density 
of the incoherent atoms at the soliton notch becomes [ac- 
cording to Eq. (30)] comparable to the condensate den- 
sity after to ~ 10 ms, see Fig. 2. 



VIII. DISCUSSION AND CONCLUSION 

It is worthwhile to compare solitons to vortices. There 
are toy model calculations in Ref. [12] for small N that 
show some incoherent atoms in the vortex core. The 
mechanism is similar to solitons: 2 atoms in the state 
with angular momentum I = 1 are depleted to states 
I = and 1 = 2. 1 = gives nonzero density in the 
core. Depletion to an empty core/notch is not limited to 
solitons. 

In contrast to solitons created by phase imprinting, 
ENS vortices are created by gradual evaporative cooling 
to a ground state of a rotating trap [13]. This ground 
state has no negative frequency anomalous mode [10]. 
This vortex ground state is stationary and has stationary 
depletion in the core. 

In conclusion, we analyze greying of a dark soliton 
due to depletion in the anomalous mode of the Bogoli- 
ubov theory. We give approximate analytic expressions 
for u{x) and v(x) of the anomalous mode valid in the 
Thomas-Fermi regime. These expressions are used to 
work out formulas for a time-dependent number of atoms 
depleted to the notch. We argue that, contrary to the as- 
sumption in Ref. [8] , the vacuum of the anomalous mode 
is far from an optimal approximation to a condensate 
with a dark soliton. There is a much better state which 
minimizes depletion from the condensate. We argue that 



this minimal depletion state can serve as a rough approx- 
imation to the state of the system right after creation of 
the soliton by the phase imprinting method. The mini- 
mal depletion state is not stationary and the density of 
incoherent atoms in the notch grows from its negligible 
initial value until it fills up the notch. For realistic pa- 
rameters like in the Hannover experiment [5] the notch 
can be filled up after just 10ms. 
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